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ABSTRACT 


Spectroscopic analysis of the cathode jet of a model coaxial magneto-plasma dynamic 
(MPD) thruster is conducted to determine electron density and temperature downstream from 
the cathode. H, line profiles were scanned from an argon-hydrogen plasma generated in the 
cathode test facility of the NASA Jet Propulsion Laboratory in Pasadena, CA. A computer 
program was written in IDL to determine the profile Doppler- and Stark half widths, which 
were used to determine temperature and electron density, respectively. Three sets of data 
from the cathode test facility were taken, while varying operating voltage, current, 
hydrogen/argon ratio, and pressure. Radial profiles for electron density and temperature were 
determined within the cathode jet. Generated plasmas ranged in electron density and 
temperature from approximately N, = 2x10'* cm? at 5000 K (0.43 eV) to 4x10'* cm” at 
15600 K (1.3 eV) . It was determined that radial density and temperature distribution within 


the cathode jet are essentially uniform. 
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I. INTRODUCTION 


A. THESIS OBJECTIVES 

The objective of this thesis is to develop a method for estimating electron 
densities and temperatures in the plasma jet of a model magnetoplasma-dynamic thruster, 
such as those presently in development for space propulsion, by analyzing the broadening 
of spectral emissions within the plasma. Both Doppler and Stark broadening effects are 
analyzed to respectively determine a radial temperature and electron density profile 
within the plasma arc. A computer program was written to analyze the broadening of 
hydrogen line profiles by unfolding and examining the Stark and Doppler half widths of 
the line profile. Three sets of data taken at the NASA Jet Propulsion Laboratory (JPL) 
Cathode Test Facility (CTF) are analyzed. 
B. THESIS OUTLINE 

This thesis is divided into six chapters and three appendices, describing the 
background, conduct, data analysis and conclusions of the experiment. Chapter I] 
provides background information for the thesis, giving a brief overview of the basic 
concepts of plasmas, magnetoplasma-dynamic thrusters, and the theory of atomic 
radiation. Chapter III continues the background discussion, outlining Doppler 
broadening, the Stark effect and Stark broadening, and line broadening analysis of 
spectral lines. Chapter IV describes the experimental facility at NASA JPL, as well as the 
data collection process. Chapter V describes the computer program written to analyze the 
data gathered, and presents and discusses the analysis of the three data sets. Chapter VI 
presents conclusions and discussion for the furthering and improvement of this and like 
experiments. 

Appendix A is the listing of the computer program, written in Interactive Data 
Language (IDL), which was written to analyze the data collected. Appendix B presents 
the results of the experiment and analysis in tabular form. Appendix C contains data 


tables used for the interpolations in the IDL program. 





Il. BACKGROUND 


A. ELECTRIC PROPULSION DEFINITION 

Electric propulsion may be defined as the acceleration of gases for propulsion by 
electrical heating and/ or by electric and magnetic body forces. As is evident in the 
definition, electric propulsion may be subdivided into three distinct categories: (Jahn, 
1968) 


1. Electrothermal. Propellant is electrically heated and thermodynamically 
expanded and accelerated to supersonic speeds through a nozzle, as ina 
chemical rocket. 


2. Electrostatic. Propellant is accelerated by the interaction of electrostatic 
fields on charged particles, such as ions or colloids. 


3. Electromagnetic. Acceleration is achieved through the interaction of electric 
and magnetic fields on a highly ionized propellant gas. 
a Plasma Thrusters 
Specifically, electromagnetic propulsion deals with the acceleration of a plasma, a 
highly ionized quasi-neutral gas, through the interaction of currents driven through the 
gas with magnetic fields established either by those currents or by external means (Jahn, 
1968). Numerous ways exist of classifying electromagnetic thrusters, from thrust mode 
to working fluid path geometry, and the names for these various thrusters are equally as 
varied. The work performed for this thesis applies to a developmental facility for a 
steady-state, coaxial, magnetoplasma-dynamic (MPD) thruster. 
a. The Coaxial MPD Thruster 
A basic diagram of a coaxial MPD thruster, showing its operation, is 
shown in Figure 2.1 on the following page. A propellant gas, commonly hydrogen, argon 
or xenon, is fed into the cylindrical thrust chamber and is ionized. Current is driven 
through the plasma by an external power supply connected across the concentric anode 


and cathode. The induction of a magnetic field by the axial current component, or 


the application of a magnetic field externally, accelerates the propellant by the resulting 
(jxB) force. 

Several models of plasma thrusters, both pulsed and steady state, have 
been flown on Satellites by the United States and the former USSR for station keeping 
and orbit mainetenance purposes (Sutton, 1992). However, MPD thrusters may have 
benefits reaching far beyond their limited application to date. Their high exhaust 
velocities and specific impulse make them appealing candidates for long duration, high 


impulse space missions, such as planetary transfers and deep space probes. 
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Figure 2.1. MPD Thruster Operation 

While some models of plasma thrusters have been flown, electromagnetic 
propulsion 1s still a largely developmental discipline. MPD thrusters employ both 
electromagnetic and plasma phenomena, both of which can prove very difficult to 
analytically model. Additionally, current designs are relatively inefficient, and large 


plasma thrusters for long duration missions have very high power requirements. (Sutton, 


122) 


B. PLASMA DEFINITION AND PARAMETERS 

Plasma, the most abundant state of matter in the universe, can be simply described 
as an electrified gas with atoms dissociated into positive (ions) and negative (electrons) 
charged particles. The degree of ionization for such a gas at thermal 
equilibrium is determined by the Saha equation: 


-UJkT 
pod 


= 2.4x10152— (2.1) 


Here, n; and n, are, respectively, the number density (per cubic centimeter) of ionized and 
neutral atoms, 7 is the gas temperature in K, k is Boltzmann’s constant, and U, is the 
ionization energy of the gas in question. 

Although all gases have some degree of ionization, not all gases may be referred 
to as plasmas. A more specific definition of a plasma is a quasineutral gas of charged 
and neutral particles that exhibits collective behavior (Chen, 1974). 

1. Plasma Temperature 

At thermal equilibrium, a gas has a distribution of velocities and, therefore, 
kinetic energies. Particle velocities most often obey a Maxwellian or Gaussian 


distribution. In one dimension, the velocity distribution is written as: 


fluydu = “e BS hy 


7 (2.2) 


where 


_ | 2kT (235) 
an 


Here, f(u)du is the number of particles (per cubic centimeter) with velocities between u 
and u+du, and v,, is the thermal velocity of the particle, dependent upon the temperature, 
T, and mass, m of the particle. 


The average kinetic energy of the particles in the gas 1s: 


ae 
oe 2 kT (2.4) 
or, in the three dimensional case: 
tieett 
a ~ 2: kT (2.5) 


Equations 2.4 and 2.5 show the most meaningful representation of temperature as 
it applies to plasmas, as a measure of average energy. Although temperature may be 
expressed in Kelvin, when dealing with plasmas it is common to express temperature in 
terms of energy, or electron volts (eV). Specifically, a 1 eV plasma means kT = 1.6x10"” 
Joules = 1 eV. Therefore, a 1 eV plasma has an E,, = 1.5 eV in three dimensions. 

an Debye Length and Sheaths 

The requirement for quasineutrality in a plasma implies that the bulk of the 
plasma be kept free of large electric potentials or fields. This is accomplished through a 
process known as Debye shielding, by which plasmas “shield out” applied electric 
potentials. 

The presence of a large potential inside a plasma would almost immediately cause 
a cloud of oppositely charged particles to surround the potential surface. This cloud of 
charged particles is known as a sheath, and it is the presence of this sheath that shields out 
the large potential that caused its formation. In a hypothetical plasma with no thermal 
motion (kT = 0) the total charge in the sheath would equal the total charge in the 
potential, and no potential would exist outside the sheath. However, in a plasma of 
temperature T, shielding will be imperfect. At the edge of the sheath, the energy of the 


electric potential is approximately equal to the thermal energy kT of the particles, which 


could escape from the potential well, allowing potentials on the order of kT/e to exist 
within the plasma. 


A measure of the thickness of this sheath, the Debye length, is defined as 


kT 
a - : (2.6) 
Anne? 





and is the distance at which the shielded potential drops by a factor of 1/e. Accordingly, 
in a collisionless plasma, the field strength as a function of distance from the potential is 


expressed as: 


d oS , mae (2.7) 


Note that 7, is used in defining the Debye length, as the electrons in the plasma, being in 
general much more mobile than the ions, do most of the shielding by creating either a 
surplus or deficit of negative charge. 

The term quasineutrality can now be defined. In order to be considered quasi- 
neutral, the effects of local concentrations of charge or the introduction of external 
potentials into the system must be shielded out in a distance short compared to the overall 
dimensions of the gas. Doing so leaves the bulk of the plasma free of large electric fields 
and potentials. Therefore, in order for a gas to be considered a plasma, Ap << L, where L 
is a measure of the size of the plasma. Furthermore, quasineutrality implies that one can 
taken, = n,=n. That is, n, and n, are considered equal, and acommon density, n, called 


the plasma density, is used. 


3: Plasma Parameter 
In order for Debye shielding to be effective there must be a sufficient number of 
particles inside the sheath. The number of particles in a “Debye sphere” of radius 4, can 


be expressed as: 
Ny = =m (2.8) 
For an ionized gas to be considered a plasma, 
Nae (2.9) 


4. Collision Frequency 
In order for an ionized gas to behave as a plasma, the motion of charged particles 
within the plasma must be controlled by electromagnetic, rather than hydrodynamic, 


forces. All plasmas have a fundamental plasma frequency, expressed as 


_ | 4nne? 
GW a —— 
: m (2.10) 


which describes the oscillation frequency of electrons inside the plasma. If t, is the mean 


time between collisions of charged particles with neutral atoms, it is required that 
wt, > 1 (2.11) 


in order for an ionized gas to behave as a plasma rather than a neutral gas. 
G: MPD PLASMA ARC 

Generating and sustaining the plasma arc in an MPD thruster entails the ionizing 
of the propellant gas as well as establishing the current and magnetic field necessary to 


accelerate the plasma jet to velocities high enough to generate sufficient thrust. The 


propellant gas, when injected into the thruster, is immediately exposed to an ionizing 
potential, existing between the cathode and the anode. The plasma is then accelerated 
past the nozzle-shaped anode. 

ie Plasma Arc Description 

In steady state operation, the plasma is readily recognizable as an intensely 
luminous column of gas extending from the cathode tip, called the cathode jet. Viewed 
externally, the exhaust plume may have a coaxial structure, consisting of the cathode jet 
and a less intense coaxial shell extending from the lip of the anode orifice (anode jet), 
separated by an even less intense region. The whole plume usually extends many anode 
diameters downstream (Jahn, 1968). Figure 2.2 shows a photograph of an MPD arc in 
vacuum tank operation. The cathode and anode jets are easily distinguishable exiting the 


anode orifice. 
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Figure 2.2. MPD Thruster Exhaust Plume (Jahn, 1968) 


D. ATOMIC STRUCTURE AND RADIATION 

The light emitted by the plasma arc described above, and indeed in all plasmas 
subjected to electrical discharges, is limited to certain discrete wavelengths that are 
characteristic of the gases that compose the plasma. Light emission from an atom occurs 


when an electron within the atom changes, or “transitions” from an initial energy level to 


a lower energy level. The wavelength of the emitted photon is 


Gal2) 


» 
tl 
<|o 


where c is the speed of light and v is the frequency of the emitted photon, which is 


directly related to the change in energy of the electron, as such: 
v= —— (2.13) 


Here, E, and E, are, respectively, the initial and final energy levels of the electron, and h 
is Planck’s constant. Furthermore, numerous models formulated to describe atomic 
structure have shown that the energy levels for electrons 1n an atom are quantized. 

1. The Bohr Model for Hydrogen 

The simplest model of atomic structure, the Bohr model, could only accurately 
describe the radiation of single electron atoms such as hydrogen or singly ionized helium. 
Nonetheless, the Bohr model accurately predicts the wavelengths of all five known series 
of the hydrogen spectrum, and works equally well on singly ionized helium. Although 
certain postulates upon which the Bohr model is based have since been disproven, the 
model serves as a suitable “jumping off point” for describing the hydrogen spectrum, 
which is of particular importance in this thesis. 

Bohr showed that the total energy of an electron “orbiting” the nucleus of an atom 


of atomic number Z could only take on discrete values 


Z 
5 = Ze ae (2.14) 
(4m€,)°2h? n? 


Here, m is the electron mass, e is the electron charge, and 1/47€, is the Coulomb’s Law 


10 


constant. h is Planck’s constant over 27m. For hydrogen, Z = 1, and 





E, = - 5 eV i eS oe (2.15) 


The integer n is the quantum number, and denotes the discrete energy level of the 
electron. n = 1, the lowest possible energy level, is considered the ground state, while 
higher energy levels are considered excited states. 

a. Quantum Numbers 

The definitions and descriptions put forth here are by no means complete. An in 
depth treatment of quantum numbers and degeneracy can be found in Eisberg and 
Resnick (1985). 

The time-independent Schrodinger equation relates the quantum state of an 
electron to its total energy: 


h? _, e 
ru W(r,8,0) + Virp(7,8,6) = Ep(r,8,0) (2.16) 


Here, yt represents the reduced mass of the system, h is Planck’s constant, and V(r) is the 
Coulomb potential of the hydrogen nucleus. £ is the total energy of the electron, as 
determined by Bohr for the quantum state n, and is called the eigenvalue. 


The terms 


gis V7 + V(r) 
2u (2.17) 


comprise an operator known as the Hamiltonian, which is designated H, such that 
Hwv(r,0,0) = Ew(r,0,0) (2.18) 


(r,8,) is the solution of the equation, called the eigenfunction. It is the product 


of the functions of the three independent variables r, 0, and , which describe the 
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position of the electron with respect to the nucleus in a spherical coordinate system: 


Wim (19,9) = RM in DP) (2. 19) 


As is evident from (2.19), the eigenfunctions are characterized by no less than three 
quantum numbers, designated n, /, and m,. Possible values for the three quantum 


numbers are as follows: 


N=). Zoe 
[ = 0,1,2,...,.n-1 
m, = -L,-leleOe Slee 


(2.20) 


The first quantum number, n, is the quantum number from the Bohr model, and 
describes the total energy of the electron. E, 1s called the eigenvalue, and n is called the 
principal quantum number. 

a. Orbital Angular Momentum 

Each electron orbit associated with a quantum number 7 has one or more 
associated values of angular momentum. The total angular momentum of the system, L, 
is quantized according to the quantum number /, termed the azimuthal quantum number, 


and 1s expressed as: 


L = yid+)n (2.21) 


This only refers to the total angular momentum of an electron orbiting a nucleus. L,, the z 
component of orbital angular momentum, is quantized according to m, , and may assume 


the following values: 


Ls = 


(Q22) 


Wz 


Orbital angular momentum of the electron is commonly represented in a 
vector model, which displays, for a single value of /, all possible states of orbital angular 
momentum as defined by m,. Figure 2.3 below shows the vector model of orbital angular 
momentum for the state / = 2. Such vector models depict the possible angular momentum 


states as vectors of magnitude L, which may be located at any point about the z axis 


defined by a cone of height +L, . 
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Figure 2.3. Vector Model of Orbital Angular Momentum for the State /= 2. Numbers are 
in units of h. (From Eisberg and Resnick, 1985) 


13 


b. Spin Angular Momentum 

In addition to orbital angular momentum, an electron posseses an intrinsic 
angular momentum, S, called its spin. Spin is quantized according to the quantum 
numbers s and ms, which may assume the follwoing values: 


Soe 
m_ = +1/2 (2.23) 


Spin angular momentum 1s quantized similar to orbital angular momentum, specifically: 


= ys(seln (2.24) 


c Total Angular Momentum 

The orbital and spin angular momenta of an electron are coupled through 
what is known as the spin-orbit interaction. This arises from an internal magnetic field 
created by the orbiting of the electron about the nucleus. This field is oriented in the 
same direction as L, and interacts with the electron’s spin magnetic dipole moment, 
which is oriented 1n the same direction as S, and produces a torque. This torque serves to 
couple the orbital and spin angular momenta, causing them both to precess about their 


vector sum, 


J=L+5 (2.25) 


Here, J is defined as the total angular momentum of the electron, and is quantized 


according to the quantum numbers j and m, , according to the quantization rules 


vIG+1)h 


J — 
J, = mh (2.26) 
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where 


o l l 
j =l+—,l-— 

ee 22) 
ee ait ye mail +7 


Total angular momentum of the electron can also be represented in a vector model, much 
like the one used to represent orbital angular momentum. 

a. Degeneracy 

From (2.20) it is readily apparent that for every value of n there are several 
different possible values of / and m,. In a one electron atom, such as hydrogen, the form 
of the eigenfunction depends on the values of all three quantum numbers, but the energy 
eigenvalue depends only upon the value of the principal quantum number, n. Because of 
this, an atom can have states with completely different behavior as defined by the 
quantum numbers, but have the same total energy. This phenomenon is called 
degeneracy, and eigenfunctions corresponding to the same eigenvalue are called 
degenerate. For example, W4) and W,,,, describe different states of a one electron atom, 
yet each state has the same total energy, E,,_, as defined by the Bohr model. In reality, the 
spin-orbit interaction results in a minute difference in energy between states with the 
same / but different 7 quantum numbers. However, the difference is very small, on the 
order of 10° eV. (Eisberg and Resnick, 1985) 

Table 2.1 below shows, for a one electron atom, all possible values of the 

quantum numbers forn = 1, 2,3. From Table 2.1, it can be seen that for every value of / 


there are (2/ + 1) degeneracies, and for every value of n there are n* degeneracies. 


; 


















Number of 
degenerate 
eigenfunctions 
for each | 








Number of 
degenerate 

eigenfunctions 
for each n 


Table 2.1. Possible Values of / and m, for n = 1,2,3 (After Eisberg and Resnick, 1985) 


3. Spectral Radiation 

An atomic electron will always seek to occupy a state with the lowest possible 
energy, within the limitations of selection rules and the Pauli exclusion principle. For the 
hydrogen atom, this energy level is the ground state, defined by the eigenfunction Wo. 
An atom in an excited state will, over time, make a series of transitions to lower energy 
levels until the ground state is reached. 

Each transition to a lower energy level is accompanied by the emission of a 


photon of energy E, - E, , and frequency 





(2-78) 


The wavelength associated with the emitted frequency is called a spectral line or line, and 
the collection of lines resulting from all possible transitions within an atom is called the 
emission spectrum of the atom. 

a. Hydrogen Spectral Lines 

The emission spectrum for hydrogen can be easily and accurately 


determined using the Balmer formula, which expresses the reciprocal wavelength of each 
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line as a function of initial and final energy states: 


K = R,(1/n, zs I/n;) (2.29) 


i 
A 


where R,, is the Rydberg constant and is defined as: 


<i acn) 7) ose = 1.097x10’m~! 
(4ne,)*4mthrc 2mhe (2.30) 





b. The Balmer Series 

The Balmer series is the family of lines arising from transitions to n,= 2 
from higher states (n = 3,4,5,...,1 = ©). The Balmer series has an infinite number of lines, 
ranging in wavelength from 6562.8 A in the visible to the series limit (1 = ©) 3645.6A in 
the near ultraviolet. 1A = 1 angstrom = 107° meters. Table 2.2 shows some of the more 
prominent lines of the Balmer series. The Balmer series is one of five known series of 
hydrogen lines, each characterized by transitions to a particular n,. The other four series 


are the Lyman (n,= 1), Paschen (n, = 3), Brackett (n,= 4) and Pfund (n,= 5) series. 





co (Series Limit) 3645.6 


Table 2.2. The a Series 


The Balmer series is of particular importance here because of its 


widespread application in the study of Stark broadening. Particular attention will be paid 


wi 


to the second line in the series, the Hg line at 4861.33 A, as it is the most widely studied 
wavelength in the theory of Stark broadening (Huddlestone and Leonard, 1965). 

C Selection Rules 

Although the Balmer formula and the Bohr model accurately predict the 
lines of all five hydrogen series, they are nonetheless incomplete models. While they 
accurately predict the change in energy and wavelength of the resulting spectral line, they 
provide very little insight into the actual atomic transition that is taking place. Modern 
quantum theory has shown that a transition can only take place if it satisfies the following 
selection rules: 


Al = 


Aj (2.38) 


1 
ou 
ie 


In order for a transition to be possible between two energy levels n,; and n, , their 
respective / quantum numbers may only differ by one, and their j quantum numbers can 
only differ by zero or one. Since the / degeneracy of the hydrogen atom causes total 
energy to be only dependent upon the principal quantum number n, the Bohr model need 
only address changes in the principal quantum number to successfully predict the 
spectrum of a one-electron atom. 

However, if the / degeneracy is removed, as it is in the presence of an 
external electric field, the total energy of a quantum state becomes dependent upon /, and 
the Bohr model and Balmer formula can no longer accurately predict the resulting 
spectral lines. The removal of the / degeneracy due to an applied electric field is the 
mechanism that causes the Stark effect, after which the phenomenon of Stark broadening 
takes its name. 

Figures 2.4 and 2.5 show two different representations of energy level 
diagrams for hydrogen, and the transitions that can result due to the applicable theories. 
Figure 2.4 shows the energy levels and transitions that are predicted by the Bohr model, 
as well as the locations of the corresponding spectral lines. Note the nonlinearity of the 


wavelength and frequency scale. Transitions that are predicted by the Bohr model but do 
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not occur are denoted by dashed lines. Figure 5 shows some of the allowed transitions in 


the Balmer (n,= 2) series, taking into account selection rules. 
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Figure 2.4. Energy Level and Transition Diagram for Hydrogen According to the Bohr 
Model. (Eisberg and Resnick, 1985) 
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Figure 2.5. Some Allowed Transitions for the Balmer Series. 
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Il. LINE BROADENING IN PLASMAS 


The two major classes of line broadening are Doppler broadening, and collision 
broadening. While Doppler broadening relates to the velocity, and ultimately, the 
temperature of the radiating atom, collision broadening is the result of collisions on the 
atomic level with other particles, either charged or neutral. The two dominant causes of 
line broadening in plasmas are Stark broadening, due to the perturbation of the energy 
states of the emitting atom by the electric fields of nearby charged particles, and Doppler 
broadening, caused by the thermal motion of the emitter. 

A. TERMS AND CONVENTIONS 

1. Line Shape 

Ideally, a spectral line is represented as a line of zero width at the wavelength at 
which it is emitted. However, even when completely isolated from its neighbors and 
shielded from all external fields, a transitioning atom will emit a line of finite width 
(Marr, 1968). The natural width of these lines tends to be on the order of hundredths or 
thousandths of an angstrom, although this is rarely the observed case. Typically, a line is 
broadened by one or more mechanisms, and takes on one of several distinct shapes 
known as a line shape. 

Line shape, or line profile, is a term used to describe the distribution of intensity 
of a spectral line over wavelength or frequency. The shape is denoted (AA) or [(Av), 


where / is the intensity distribution. 


a 


Ne va, (3.1) 


i, and v, are the wavelength and frequency of the emission, also called the center 


wavelength and frequency. /(AA) can take on one of three different shapes. 


Py 


a. Lorentzian Profile 
Also known as a “dispersion-type” (Huddlestone and Leonard, 1965) 
shape, the Lorentzian shape is the natural line shape of an unbroadened spectral line 


(Marr, 1968), and is represented by the function 


| 
Ai? 
AX 


(AX) « 
] + 





(3.2) 


Ak, ,, the half width of the profile, is the wavelength difference at which the intensity 
drops by a factor of two. 2AA,, is known as the full half-width of the profile. Lorentzian 
profiles are associated with collision broadened lines, of which Stark broadening is a 


type. A lorentzian profile is shown below in Figure 3.1. 
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Figure 3.1. Lorentzian Profile, Ax,, = 0.707. 
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b. Gaussian Profile 
A Gaussian shape, also known as a Maxwellian distribution, takes the 
form 


{44} 
KAA) = I,e (3.3) 


where J, is the intensity at the line center, and 1, the Doppler width, is the wavelength 
difference at which the intensity is reduced by a factor of 1/e. The half width of a 
Gaussian line shape is 


(3.4) 
Ad.) = Adpyin2 


The term Doppler width is used here in reference to Doppler broadening, which produces 
line shapes that are distinctly Gaussian (Griem, 1964). Doppler broadening will be 
discussed in detail later in this chapter. An example of a Gaussian profile is shown in 


peure 3.2) 


y = exp[-(x/2)*2) 





Figure 3.2. Gaussian Profile 
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C. Voigt Profiles 

In many instances, a spectral line will exhibit neither a purely Lorentzian 
nor a purely Gaussian shape, but will instead show features of both distributions. Such a 
profile is called a Voigt profile, and is formed by the convolution of a Gaussian with a 


Lorentzian profile: 


WMA) = [-1,(Ad)1,(AA*-AA)d(AA) (3.5) 


—o 


The resulting profile is expressed as a function of the dimensionless variable w (Marr, 


1968): 


2 





Ad) « [f° __&@" sg 
i az + (w-vY ‘ oe 
where 
: mee 
os Ai, (3.7) 


and a is called the damping ratio and is defined as the ratio of the Lorentzian to the 
Doppler half widths: 
- Ai, 

Ad, vin2 (ue) 


a 


A Voigt profile will be the resultant line shape when both collision-type and Doppler 
broadening are present in a line profile. Examples of Voigt profiles with varying 


parameters can be seen in Figure 3.3. 
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Figure 3.3. Voigt Profiles for Damping Parameter (a) of 0, 0.5 
B. DOPPLER BROADENING 

Doppler broadening takes its name from the so-called Doppler effect, which states 
that radiation emitted from a source moving towards or away from an observer 1s shifted 


in frequency by 
ice) 


where Vv, is the frequency of the source at rest, v, is the velocity of the source, and c 1s the 
speed of light. Radiation from a source moving toward an observer experiences a 
positive shift, while a negative shift occurs when the source is moving away from the 
observer. The resulting wavelength shift is also expressed as 


VY 
aaa 
dunia hey S (3.10) 
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where A, is the unshifted wavelength. In a plasma, the random motions of radiating 
particles causes a Doppler shift in the emitted radiation, which manifests itself in a 
phenomenon known as Doppler broadening. Assuming the motion of particles within the 
plasma is of a thermal nature, the velocity distribution of the emitting atoms is Gaussian, 


and 


== (Gal 


represents the fraction of particles moving in the line of sight with velocities between v, 


and v, + dv, , and v = (2kT/m)” being the most probable velocity. By substituting AA for 
v, from (3.10), and defining 





v 
Ai, = ri (3.12) 
(3.11) may be rewritten as: 
d 1 {=} 
un Ee (3.13) 
N Avr 


Lastly, assuming the total intensity, /, is proportional to N, and (AA)d(AA) is proportional 
to dN, the formula 


{arp 
(AX) = —_ i 


dV % 


(3.14) 


is arrived at for the line shape of a Doppler broadened line. 
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It is apparent from 


~ 


o | 2k 


Cc m 


hn] 


Ai, = 


D Bs15) 
that if the Doppler width of a profile is known, the plasma temperature, at least for the 
emitting atoms, can be determined. Doppler broadening can also result from the bulk 
motion of the plasma, and is evident in many pulsed plasmas (Huddlestone and Leonard, 
1965). 
C. STARK BROADENING 

Unlike Doppler broadening, Stark broadening is only a single category of what is 
known as collision or pressure broadening. Collision broadening is largely due to the 
interaction of the emitter with surrounding particles, and is subdivided into three distinct 
types: 

1. Resonance broadening due to interaction with atoms of the same kind. 


2. Van der Waals broadening due to interaction with neutral molecules or atoms 
of different species. 


3. Stark broadening due to interaction with charged particles (ions and electrons). 


Named after the so-called Stark effect, Stark broadening is the broadening of the natural 
(Lorentzian) shape of a spectral line due to the interaction of the emitter with the 
Coulomb forces of nearby ions and electrons, which are present in abundance in dense 
plasmas. 

1. Stark Effect 

The Stark effect takes its name from the investigations of Stark in 1913 on the 
emission spectrum of hydrogen in the presence of an electric field (Marr, 1968). Stark 
discovered that a splitting of spectral lines occurs when the radiators are in the presence 
of an electric field, due to the removal of degeneracies in the unperturbed eigenvalues of 
an atomic electron. The quantum mechanical treatment of the Stark effect 1s rather 


complex compared to other splitting mechanisms such as the Zeeman effect, and is not 


yee 


needed in the analysis of a spectrum (Herzberg, 1944). The theory will therefore be 
presented in its elementary form, with the appropriate references made to quantum 
mechanical conventions. 

a. Elementary Stark Effect Theory 

In the presence of an electric field of strength F, the magnetic moment 
associated with the total angular momentum, J, of an atom is unaffected. However, the 
field will polarize the atom, by displacing the positively charged nucleus from the center 
of gravity of the negative charge distribution. This results in the creation of an electric 
dipole moment proportional to F, say aF, where a depends upon the orientation of the 
electron charge distribution to the field, and is therefore a function of J. This is shown in 
Figure 3.4 where the nucleus N is shown separated from the center of negative charge S$ 
(Marr, 1968). 

A precession of J about the electric field direction occurs due to the 
interaction of the electric dipole with the electric field, such that the component M, in the 
direction of the field remains constant. The rate of precession increases with increasing 
F, and the energy change for given angular momentum state J is given by the product of 


the field strength and the electric dipole, 


- 2 
AE = KM,F 3.16) 


where K is a constant denoting a fractional portion of M, (Herzberg, 1944). 

It can be seen from (3.16) that the energy shift due to the Stark effect is 
proportional to the square of the field strength. This is termed the quadratic Stark effect. 
Also, it is of note that term components that differ in the sign of M, have the same energy, 
as reversing the rotation of an orbiting electron (changing M, to -M, ) does not alter the 


electric dipole moment. 
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Figure 3.4. Vector Diagram for Electric Field F Coupling for Atoms. S is the center of 
gravity of the electrons with resultant angular momentum J perpendicular to the orbital 
plane. (After Marr, 1968) 


b. Stark Effect for Hydrogen 
Unlike other atoms, hydrogen exhibits a strong Jinear Stark effect, that 1s, 
the energy shift increases linearly with field strength. If a perturbing electric field of 


Strength F is applied to an atom, the potential of that field is: 


eFz = eFrcosO (3.17) 


Time independent perturbation theory states that the total energy of a state n under the 


influence of a perturbing potential is 


liege |e 
ea ea ya ne, (3.18) 
" "(B, ~ E,) 
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where H’,,, is the interaction Hamiltonian between energy states n and n’, defined as 


ee - [Ye v¥ AY (3.19) 


(3.19) defines the interaction between the two states n and n’ due to a perturbing 
potential, v. In most atoms subject to electric fields, the first order term (H7’, ,) of the 
perturbation vanishes, and the quadratic term must be used to a first approximation. 
However, if states are degenerate to the extent that two states with different / quantum 
numbers have the same energy (as in the hydrogen atom) it is possible to have a nonzero 
first order effect provided that AJ = +1 and Am # 0 (Marr, 1968). 

Both the Bohr theory and modern quantum mechanics show that, for 
hydrogen, an energy level n splits into 2n - 1 equidistant energy levels, the energy 
difference of which increase linearly with field strength (Herzberg, 1944). 

pe Stark Broadening Theory 

In the presence of a significant quantity of ions and electrons, as in a plasma, a 
radiator is subjected to electric fields which consequently perturb its normal energy levels 
by way of the Stark effect. Varying electric field strengths within the plasma cause a 
spreading of the Stark levels in individual atoms, resulting in a statistical broadening of 
the spectral lines, known as Stark broadening. In dense plasmas, Stark broadening is the 
predominant form of line broadening. 

If the perturbing ion or electron is treated as a classical particle, then the duration 


of each perturbation in a statistical distribution of perturbing particles is defined by 

t= PI (3.20) 
where £ is the impact parameter of the perturbing particle, and v,, is the average velocity 
of the perturber. Due to the large difference in the velocities of the fast moving electrons 


and slow ions, t, varies greatly for the particles. This resulted, initially, in two theories of 


Stark Broadening, called the impact approximation and the quasi-static approximation. 
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a. Quasi-Static Approximation 

The quasi-static approximation 1s used in calculating the Stark broadening 
due to ions, as T, 1s generally larger than the average time between collisions with ions. 
For that reason, the motion of ions may be neglected completely, and a static Stark effect 
is assumed, taking into account the statistical distribution of electric fields at the radiating 
atom (Marr, 1968). 

In computing a line profile using the quasi-static approximation, the 


splitting of the line resulting from the Stark effect is first calculated, which is defined as: 


x u 
Vie Cit (3.21) 
where p = | for the linear Stark effect, 2 for quadratic, and Ci, is the Stark coefficient for 
transitions between initial state 7 and final state fin a radiating atom, and is derived from 
applying perturbation theory (Griem, 1964). The line splitting is then averaged over the 


probability distribution of all different electric microfields within the plasma, denoted 


beatae 
WwW = eee —L— = ————— 
HF) F A e F y(B) (3.22) 


F,, called the normal or Holtsmark field strength (Griem, 1964), depends only upon ion 
density and is defined as | 


ee SIO N22 Vic 


° Ane, (3.23) 


N 1s the ion density of the plasma. 

b. Impact Approximation 

The impact approximation for Stark broadening must take into account the 
infrequent and short-lived interactions of the radiator with highly mobile electrons, and 


therefore cannot assume a time-independent perturbation as in the quasi-static 
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approximation. Determining an impact-broadened line profile first requires the solution 


of a time-dependent Schrodinger equation of the form 


mM = [H, + ViOlx) (3.24) 


where H, is the Hamiltonian of the unperturbed radiator, and V,(Z) is the potential of the 


perturbing electron. x(t) describes the state of the emitting electron, and is defined as 
X(t) = T,(t,0)x(O) (3.25) 


where 7,,(t,0) is called the time evolution operator, which relates wave functions at time f 
to those at t= 0 (Griem, 1964). 

This equation must be solved for 7, for all V,(t) produced by perturbers passing 
with various velocities, impact parameters, and times of closest approach. These 
solutions are then substituted into a correlation function which 1s statistically averaged 
over all possible perturber trajectories. Finally, the Fourier transform of the correlation 
function is taken in order to obtain the impact broadened line profile. (Griem, 1964) 

a: Stark Profiles 

In order to generate a Stark broadened line profile using both the impact and 
quasi-static approximations, the quasi-static theory is first applied, subjecting the 
radiating atom to a Static field causing energy level shifts and line splitting. After that, 
the electron broadening is superimposed upon the resulting Stark patterns (Huddlestone 
and Leonard, 1965). Stark profiles have been generated and tabulated, most notably for 
the Balmer series of hydrogen, for use in plasma spectroscopy. Stark broadening 


produces a Lorentzian line profile. 


Se 


The tabulated profiles themselves are one half of a symmetrical profile, and are 
generally presented as functions of S(a) versus «. a is known as the reduced wavelength, 


and is expressed as 





Zz 
et cubes (3.26) 
F 2ncF 


F) is the Holtsmark field strength, dependent upon ion density. From S(q), a profile in 


the wavelength scale can be obtained from 


: da _ p-l 
I(A) = S(a) a Fy, S(a) (3.27) 


Stark profiles for hydrogen have been extensively calculated using numerous 
approximations. The most accurate Stark tables for hydrogen have been generated using 
the “Unified Classical Path Theory”, often called the “Unified Theory” of Stark 
broadening, developed by Vidal, Cooper, and Smith (1970). While adhering to the quasi- 
Static approximation for ions, the Unified Theory of Stark broadening develops a more 
complex model for electron collisions. This gives better agreement with experimental 
data than the impact theory does, particularly in the line wings and transition region 
between the line wings and line center. Stark profiles for a given spectral line are 
normally presented in tabular format according to electron density and temperature of the 
plasma environment of the radiator. The profiles are well known for several lines of the 
Balmer series, and have been computed for other hydrogen series as well, most notably 


the Lyman series. All computed profiles are normalized such that 


[sada = 1 (3.28) 
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The most notable and distinguishing aspect of Stark broadened line shapes occurs 
at relatively high electron (ion) densities (NV, 2 3 x 10'*cm”’) where the Stark effect due to 
ions is great enough that the spectral line splits. When this happens, the line appears as 
two Lorentzian lines side by side, which overlap at the line center. Figure 3.5 shows the 
comparison of theoretical and experimental Stark profiles for the Hg line at an electron 
density of 6.4x10'® cm®. Generally, the splitting of the line is stronger in theory than in 
experiment (Vidal, Cooper and Smith, 1972), although in Figure 3.5 it is quite evident in 


both the measured and theoretical profiles. 
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Figure 3.5. Comparison of Experimental and Theoretical Stark Profiles of the H, Line. 
N, = 6.4x10'° cm®, T, = 12,210 K. (From Huddlestone and Leonard, 1965) 


4. Electron Densities 
The most common use of Stark broadening as a diagnostic tool is in determining 


electron densities in plasmas (Huddlestone and Leonard, 1965). The most common way 
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of doing this is by comparing the full half width of a spectral line with the half widths of 


computed Stark profiles. Electron density is then determined using the relation 


° 


Akin = 2F hip = 2.5x10° a, ,Ne- A (3.29) 

where Ad, is the full half width of the line profile, and «,, is the value of « at which 
S() is one half its maximum value. o,,. must be retrieved from the theoretical profile. 
This has already been done for various hydrogen lines, at various electron densities and 
temperatures, by Griem (1974). J, is the electron density of the plasma. Using half 
width comparison, errors on the order of 10 - 15% can be expected, the lower limit being 
achievable for lower ( 10'° cm”) electron densities (Griem, 1964). 

The second way of determining electron densities is by a least squares fit of an 
observed line profile to the theoretical Stark profiles. This method is the better of the two 
for higher electron densities (~10'° cm”) where splitting of spectral lines (especially the 
H, line) becomes significant. At lower electron densities ( 10’ cm”), error should be 
on the order of 5%, while an upper limit of 10% error can be expected for higher electron 
densities (Vidal, Cooper and Smith, 1972). 

D. ELECTRON DENSITY DETERMINATION IN PLASMAS 

Very rarely will a spectral line emitted in a plasma be either purely Doppler- or 
Stark broadened. In fact, this is almost never the case; both forms of line broadening will 
invariably contribute to the line shape, and a Voigt profile will be the result. Before Stark 
broadening analysis can be applied to the spectral line, the Voigt profile must be analyzed 
to determine the extent to which each broadening mechanism is present, and the Doppler 
broadening must be factored out of the line profile. 

Ie Half Width Estimation 

When concerned only with the half-width of the observed line profile, it is 
sufficient to determine the Voigt parameters which best fit the observed profile, from 


which the Lorentzian and Gaussian half-widths can be extracted, and density and 
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temperature can be determined. The relationships between the Stark, Doppler, and 
observed half-widths of line profiles are presented in Appendix C. 

Once the Stark half width of the profile is CRON can be compared to the 
known half widths of theoretical Stark profiles, which are tabulated as functions of N, and 
T in Appendix [J of Griem (1974) for several hydrogen and ionized helium lines. Given 
a reasonable value of 7, either measured or determined from the Doppler width of the 
profile, and the Stark width of the observed profile, an electron density can be 
interpolated from the half width tables presented in Griem (1974). Figure 3.6 shows, for 
H,, electron density versus half width for T = 5000 and 40000 K . 
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Figure 3.6. Electron Density Versus Full Half Width of the H, Line for T= 5000 and 
40000 K. 
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IV. EXPERIMENT SETUP 


A. CATHODE TEST FACILITY 

The cathode test facility (CTF) at NASA-JPL is designed to observe and analyze 
cathodes under the simulated operating conditions of an MPD thruster. It contains all the 
necessary hardware and software to observe and record line profiles emitted from the 
plasmas generated in the facility. A complete description of the CTF is presented in 
Goodfellow (1996). 

1. Vacuum Chamber 

The main component of the cathode test facility at NASA-JPL is a stainless steel 
vacuum chamber, 0.5 m in diameter and 2.4 m long, consisting of four water-cooled, 
cylindrical segments. A water-cooled, copper liner is inserted in the middle two segments 
to facilitate long duration cathode operation. 

The cathode fixture is mounted on the vacuum chamber door, and consists of two 
coaxial tubes, electrically isolated from each other and the door. The cathode itself is 
clamped to the inner tube, which serves as the current feed for cathode operation. The 
cathode protrudes through a water cooled copper disk mounted on the end of the outer 
tube. The propellant gas is injected between the two tubes, and enters the discharge 
chamber (the first segment of the vacuum chamber) through an annulus. formed around 
the base of the cathode. The cathode itself is 76 mm in length and 9.5 mm in diameter 
with a hemispherical tip, and is made of 2 percent thoriated tungsten. 

The anode is a water cooled copper ring 7.6 cm in diameter, which is mounted on 
a flange between the first two segments of the vacuum chamber, from which it is 
electrically isolated. The last segment of the vacuum chamber contains a heat exchanger 
composed of water-cooled, finned copper tubing, which cools the exhaust before it enters 
the pumping system. Several optical ports line both sides of the chamber, providing 
access to the discharge chamber (particularly the cathode region itself) and the plasma 
plume. The rear of the chamber has a port which allows viewing of the cathode tip and 


plume along the axis of the tank. 


Sif 


The chamber is pumped down using a 610 liter/sec Roots blower, backed by a 
140 liter/sec Stokes mechanical pump. Vacuums of less than 0.13 Pa (0.975 mTorr) can 
be achieved without propellant flow into the chamber, while a vacuum of approximately 
80 Pa (0.6 Torr) can be achieved with an argon flow rate of 0.75 g/s (Goodfellow, 1996). 
Typical operating pressures during the experiment ranged from 11 Torr (1.467 kPa) to 30 
Torr (4 kPa). Pressure in the discharge chamber was controlled either by adjusting 
propellant gas flow rate into the chamber, or by adjusting a gate valve between the 
chamber and the pumps. 

Power to the arc is supplied from two Miller welding power supplies, each 
capable of supplying a continuous 1500 A at 60 V, or 2000 A at 50 V for up to 20 
minutes. The initial arc breakdown is achieved with a 4 A, 850 V start supply. 

2 Tank Instrumentation 

Arc current is monitored with precision shunts replacing the factory shunts on the 
Miller power supplies. Terminal voltage is measured at the current feedthroughs into the 
vacuum tank. 

A Sierra Instruments Side-Trak Model 830 flow meter and a Micromotion Model 
D6 flow meter are used to monitor propellant flow into the thrust chamber. Propellant 
flow is controlled with a throttling gate valve located upstream of the inlet to the cathode 
fixture. An MKS Baratron capacitance manometer (0 - 133 Pa range) is used to monitor 
tank pressure. 

The above parameters, as well as various system temperatures (cathode, anode, 
etc.) were monitored and recorded with a Macintosh computer utilizing LabView 
software and Opto-22 data acquisition hardware. 

3. Data Collection Equipment 

A system for performing emission spectroscopy within the plasma plume was 
already in place in the cathode test facility. A lens was mounted outside an optical port at 
the thrust chamber, which focussed an image of the cathode on a screen. The inlet to a 
length of 100 micron diameter fiber optic cable was located at the center of the screen, 


flush with the surface. The image screen and fiber inlet were both mounted on a 
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micrometer-operated translation stage, which allowed the fiber inlet to be placed at any 
point within the image. The translation stage allowed adjustments in position to less than 
0.5 mm. 

Figure 4.1 shows a basic schematic of the data collection system. Light exiting 
the fiber was focused on the entrance slit of a one meter McPherson scanning 
monochromator. A 1200 groove/mm grating was used in the monochromator to disperse 


the light, and the variable entrance and exit slits were typically set at 30 microns. 
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Figure 4.1. CTF Data Collection Equipment Setup. 
Spectral resolution of the spectrometer is on the order of | A. The output of the 


monochromator is amplified and stored as a voltage in a Macintosh computer, again using 


LabView software. 
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B. DATA COLLECTION 

Data were gathered while the CTF was operating on an argon plasma, with small 
amounts of hydrogen (10 to 15% by volume) intermixed to allow observation of the H, 
line. Under normal operating conditions, the plasma arc flows from the tip of the 
cathode, forming a cathode jet as defined in Chapter Il. Sets of hydrogen line profiles for 
analysis were gathered from within this cathode jet. The scanned wavelength interval 
was from 4856A to 4866A, to allow inclusion of as much of the wings of the line as 
possible without serious overlap of nearby argon lines. 

Initially, the horizontal micrometer translation stage was adjusted so that the inlet 
to the fiber optic cable was located directly off the tip of the cathode in the image plane. 
From this location on the cathode axis, the image plane could be translated vertically, so 
that the fiber inlet moved radially within the plasma plume. The vertical micrometer was 
adjusted in increments of 0.083 in (2.108 mm) between readings, which were taken in 
sets of 31. One reading was taken on the cathode axis, and 15 readings were taken to 
either side, up or down, of the cathode. Thus, a set of 31 readings covered a “sliver” of 
the plasma jet 2.49 in (63.246 mm) in diameter, centered about the cathode axis. 

1. Data Format 

As stated above, a 10 A interval was scanned about the location of the Hg 
(4861.33 A) line, from 4856 to 4866 A at a resolution of 500 data points per angstrom. 
Each scan was stored on a Macintosh personal computer, using LabView software, as a 
file consisting of two arrays - wavelength, recorded in angstroms, and intensity, recorded 
as a voltage level. Each array consists of 5000 double precision floating point data 
values. For bookkeeping purposes, each set of scans was given a group name, and 
operating conditions such as argon/hydrogen flow rates, pressure, voltage and current 
were also recorded for each set of data. The cathode temperature profile was also 


monitored during data collection. 
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Figure 4.2. Scan Location With Respect to Cathode Tip. 
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V. ANALYSIS 


In order to determine electron temperatures and denisities from the scanned line 
profiles, a computer program to determine characteristic line widths (both Doppler and 
Stark) was written in Interactive Data Language (IDL). Sequentially reading all 31 files 
in a data set, the program attempts to curve fit a Voigt profile to each line in order to 
estimate the parameters needed to extract the appropriate line widths. The appropriate 
characteristics are then determined either by a simple formula in the case of temperature, 
or by interpolation of tabulated theoretical data in the case of electron densities. Three 
separate sets of data were analyzed, varying in distance from the cathode tip and 
Operating parameters, for which an overview of the analysis is given in this chapter. A 
documented listing of the IDL program is presented in Appendix A, and tabulated results 
of the analysis are presented in Appendix B. 

A. ANALYSIS PROGRAM 

The analysis program consists of eleven procedures which perform the tasks of 
reading and preparing the data, determining the contributions of Doppler and Stark effects 
to line width, calculating the electron temperatures and densities, and storing the results 
in a single data file for a set of scans. The majority of the procedures are used in the 
fitting of the data and determination of line widths. The program is initiated with a single 
call to a master procedure, called LINBROAD, which executes the above described 
functions in order. In turn, LINBROAD executes three sub-procedures: 

1. FITFILE, which fits the appropriate Voigt profile to the scan and determines 

the Stark and Doppler widths of the line profile. 


2. NETE, which takes the line widths output from FITFILE and computes the 
electron temperature and density of the profile. 


3. STORIT, which writes the results to a data file with a user-specified name. 
The first two procedures in turn call one or more sub-procedures of their own to 


accomplish their tasks. 
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1. FITFILE 

FITFILE is an iterative procedure which first reads each data file one at a time 
using a procedure called GETFILE. Following the input of a single data file, FITFILE 
successively fits Voigt profiles to the data, adjusting the guessed Doppler width while 
keeping track of the Chi-Squared parameter of the fit, until a minimum Chi-Squared is 
reached, indicating a best fit. 

Once a best-fit Voigt profile is found, the HALFWIDTH procedure is called in 
order to determine the half-width of the fit. Following that, the GETPAR procedure is 
executed to determine the fractional contribution of both Doppler and Stark broadening to 
the line width. The PROWIDTH procedure is then executed, which determines the full 
half width of the data. Lastly, FITFILE computes the Doppler and Stark widths of the 
data from the outputs of the GETPAR and PROWIDTH procedures. 

a. GETFILE 

GETFILE first reads the data file and chops off the first and last 200 data 
points in each array (wavelength and intensity), as the extreme wings of the scanned 
profile contain significant portions of neighboring argon lines which would otherwise 
adversely affect the fitting of the data. The procedure PROFCOR is called, which 
initially fits a Lorentzian profile to the data in order to determine the location of the line 
center. PROFCOR also determines an offset parameter, as the voltage (intensity) levels 
of the scan settle into negative values out in the line wings, which must be corrected for. 

After the intial fit, the profile is further reduced in size to 3000 data points, 
and the profile spans 3A to either side of the line center. The profile is then “folded” in 
half about the line center, averaged and then “unfolded” to produce a smoother, less 
“noisy” profile to accomodate the fitting of a Voigt profile. 

b. GETPAR 

Once a Voigt profile has been fitted to the data and the half width of the 
Voigt profile is determined, GETPAR is called to determine the relative contributions of 
Doppler and Stark broadening, based upon the Doppler width and half width of the fitted 


Voigt profile. This is done by executing a second-order polynomial fit of the ratio of the 


te. 


Lorentzian half width to the full Voigt half width as a function of the ratio of Doppler 
half width to the full Voigt half width. The data to be fitted is presented in Huddlestone 
and Leonard (1965), and is reproduced in Appendix C. 
e PROWIDTH 
PROWIDTH determines the full half width of the observed profile by first 
dividing the data into bins of 20 points each and taking the average of those bins, 
reducing the profile from 3000 to 150 points, and then finding the location of the bin on 
either side of the line center that most closely represents the half width. The right and left 
side half widths are then added to produce the full half width of the profile which, along 
with the ratios determined in GETPAR, is used to compute Doppler and Stark full half 
widths. 
2: NETE 
Given values for the Doppler and Stark widths of the line profile, NETE 
determines first the temperature, and then the electron density of the plasma. The 
temperature determination is easily accomplished by manipulating equations (3.4) and 


(3.15), giving 
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where AA,” is the Doppler half width. 

Electron densities are determined by interpolation of a table of theoretical Stark 
widths (in log space) as a function of temperature, as presented in Appendix III of Griem 
(1974). The portion of this table used is presented in Appendix C. The Stark widths for a 
given electron density as a function of electron temperature are fitted with a cubic spline, 
which is then interpolated for the temperature determined in (5.1). The spline is 
performed in log space with respect to temperature, resulting in the log of the electron 


density as a linear function of the log of the Stark width. A linear interpolation is then 


45 


performed for the log of the observed Stark width of the profile. The inverse logarithm of 
this interpolated value yields the electron density for the profile. 

>: STORIT 

After the electron densities and temperatures have been determined for a data set, 
STORIT writes the density and temperature vectors, as well as the full half width, 
Doppler and Stark widths, and the ratios of the Doppler and Stark widths to the full half 
width to a file with a user-specified name. STORIT also records the distance of each 
measurement from the cathode axis. 

B. DATA REPRESENTATION 

Each set of data is given a name for identification, consisting of the letters “CTF”, 
for the cathode test facility at JPL, followed by a number which designates, in sequence, 
the running of the CTF for that purpose. For example, the data set CTF144 denotes the 
144th running of the facility for gathering data. The three data sets analyzed here are 
CTF’s 144, 145, and 146. 

For each data set, operating conditions which may have some bearing on the 
analysis results are recorded. These are argon and hydrogen flow, in SLPM, tank 
pressure in Torr, and operating voltage and current. 

Scans within a data set are numbered sequentially, from 1 to 31, with the file 
name “scan1", “scan2", etc., “scan1" being located on the cathode axis. Scans 2 through 
16 make up the 15 scans above the cathode, and scans 17 through 31 lie below. 

C. DATA ANALYSIS RESULTS 

Of the three sets of data analyzed, one, CTF144, was taken at a location just off 
the tip of the cathode, while the remaining two were taken at a distance of approximately 
one half inch (1.27 cm) from the cathode tip. Tank pressure, hydrogen to total volume 
flow ratio, voltage and current were similar for CTF144 and CTF145, with CTF146 being 
taken at slightly higher values. While differing in numerical values, electron density and 
temperature behavior across the plasma plume were similar for all three data sets. Table 
5.1 summarizes the operating conditions and extreme values for electron density and 


temperature for all three data sets. 
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CTF144 CTF145 CTF146 
Distance from Ocm ee rem 1.27 cm 
Cathode (0.5 in) (0.5 in) 
Argon/Hydrogen 6.08 Ar 6.07 Ar 6.12 Ar 
Flow Rate 0.8 H, 0.9 H, 1.0 H, 

| Pressure 12.25 Torr 11.8 Torr 21.4 Torr 
a a a 
1400 A 1400 A 1400 A 

Max T 8600 K 7200 K 15600 K 

Min T 6200 K 5000 K 12900 K 


Max N, 2.63x10" cm? 2.41x10"4 cm? 3.95x10"* cm? 
Min N, 2.30x10"* cm” 2.14x10"* cm? 3.61x10'* cm? 


Table 5.1. Summary of Operating Conditions and Analysis. 














Figures 5.1 and 5.2 below show the radial electron density and temperature 
profiles of all three data sets. All results were obtained using LINBROAD. Radius from 
the cathode axis is represented on a linear scale and electron density and temperature on a 
logarithmic scale. Apart from small fluctuations in N, and T,, all three data sets are 
reasonably well behaved. They suggest that the density and temperature distributions 
across the cathode jet were nearly constant. No apparatus presently exists at the NASA- 
JPL cathode test facility for direct measurement of electron densities and temperatures 


within the vaccuum chamber, and so these results could not be checked independently. 
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x - CTF144: 12.3 Torr, 11.6% H, 17.6 V, 1400A 
o - CTF145: 11.8 Torr, 12.9% H, 17.8 V, 1400A 
*-CTF146: 21.4 Torr, 14.0% H, 22.9 V,1400A 
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Figure 5.1. Radial Electron Density Profiles Obtined for CTF’s 144, 145 and 146 Using 
the LINBROAD Program. Tank pressure, hydrogen to total volume flow rate, cathode / 
anode voltage and cathode current are listed for each data set. 
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X - CTF144: 12.3 Torr, 11.6% H, 17.6 V, 1400 A 
0 - CTF145: 11.8 Torr, 12.9% H, 17.8 V, 1400 A 


*- CTF146: 21.4 Torr, 14.0% H, 22.9 V,1400A 
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Figure 5.2. Radial Electron Temperature Profiles Obtained for CTF’s 144, 145 and 146 
Using the LINBROAD Program. Tank pressure, hydrogen to total volume flow rate, 
cathode / anode voltage and cathode current are listed for each data set. 

1G CTF144 

Data were first collected for R > 0. Then, the fiber optic inlet was centered again, 
and data was collected for R <0. During the collection of the lower half of the data set, 
operating conditions within the vacuum chamber were observed to change. An increase 
of hydrogen flow from 0.8 to 1.4 SLPM was observed, along with a major upward shift in 
the cathode temperature profile near the tip. This shift could be responsible for the slight 
change in densities and temperatures shown in the region of -1 < R<0Ocm in Figure 5.1. 
It appears as if electron densities and temperatures reached a sort of “steady state” 
condition towards the end of the recording of the data set. Figure 5.3 below shows the 


center scan (R = 0 cm) from this data set, before being analyzed. Note the two nearby 
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argon lines, one at 4866A and the other, weaker line beginning to show at around 4856A. 
It can be seen that a fair amount of noise is present in the wings of the 4861.33A | 

hydrogen line,as well as at the line center. The scaling of the intensity into both positive 
and negative voltage values is also evident. It was this offset that necessitated an “offset 


constant” be included in the curve fitting programs. 


6 FWHM = 0.85 Angstroms 


5 T= 7000'K 


4 Ne = 2.35e+14 cm4-3 
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Figure 5.3. CTF144 Center Scan, Raw Data. 


a CTF145 
On the average, the CTF145 lines were narrower than those in CTF144, despite 
similar operating voltage and current. Hydrogen flow rate was slightly higher than 


CTF144, with a difference in pressure of -0.5 Torr. No shifts in operating conditions 
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were observed during the data collection, and again temperature and density behavior is 
fairly constant throughout the profile. Maximum and minumum density values differ by 
12.6%, and density values are between 7 and 8% less than those observed in CTF144. 
Temperatures are between 16 and 19% less than those of CTF144 as well. The decreased 
operating pressure during this data collection may contribute to this “falling off’ of 
density and temperature. Being that the cathode jet is a fairly well contained plume of 
plasma, both by magnetic and tank pressure, it is doubtful that the location of this data set 
a mere one half inch from the cathode tip would show any major differences in electron 
densities of temperatures than one taken directly off the cathode tip. 

Figure 5.4 shows the center scan for CTF144. The realtive intensity of the 
accompanying argon lines is noticeably decreased. The higher hydrogen flow rate (12.9% 
of total vice 11.8%) may partially account for this by increasing the number of hydrogen 


emitters relative to argon. 
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Figure 5.4. CTF145 Center Scan, Raw Data. 
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3 CTF146 

Scan 30 of CTF146 could not be completely loaded from disk due to a disk error, 
so the data set was shrunk to 27 data files by eliminating the last two scans on either end 
of the data set. For this set, tank pressure and voltage were significantly increased, to 
21.4 Torr and 22.9 volts respectively. The resulting line profiles are significantly broader, 
reflecting much higher temperatures and electron densities compared to the other two 
data sets. Under the increased pressure the electron density distribution is much more 
uniform, with maximum and minimum values differing by only 9.4%, with no increasing 
or decreasing trends evident across either the temperature or the pressure profile. With 
increased power, (V/) the hotter cathode is stripping more electrons and generating more 
ions, while the increased pressure acting on the cathode jet forces the particles in the 
plume closer together into a denser, more uniform distribution. 

The CTF146 center scan, shown in Figure 5.5 1s very noticeably broadened, 


especially in the line wings and the “transition region’”’ between the wings and line center. 
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Figure 5.5. CTF146 Center Scan, Raw Data. 
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D. ANALYSIS COMPLICATIONS 

The scaling of the intensity values of the line profiles presented problems in the 
analyzing of the data, particularly in the fitting routines used to model the data. As is 
evident in Figures 5.3, 5.4, and 5.5, the voltages representing line intensities are scaled 
such that the line wings settle out at negative levels. Statistical distributions, by their 
nature, do not accommodate negative values, and trying to curve fit a Voigt profile to one 
of these data sets is difficult, requiring the inclusion of an “offset parameter’ to estimate 
the actual (positive) background levels with which the profile merges in the line wings. 
Attempting to fit the profiles with the LORENTZ function written into the program posed 
little problem, outside of ambiguities in background intensity levels resulting from the 
estimation of the offset. However, the VOIGT function inherent to the IDL language ran 
into problems when attempting to fit profiles stretching into negative voliages, and even 
with an offset parameter figured into the fit, the initial guesses for fit parameters had to be 
carefully chosen to keep the routine from producing meaningless output. 

Attempting to first offset the data into positive voltages before attempting to fit 
the data proved to be useless, and no good estimate of background levels could be arrived 
at outside of the offset parameter generated by curve fitting functions. Though reasonably 
accurate, the fits generated for the profiles were not as good as expected, and some small 


amount of accuracy was no doubt sacrificed. 
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VI. CONCLUSIONS 


As stated above, no means of directly measuring electron densities are in place at 
JPL, and this work can only be considered a first (and so far, only) estimate as to the 
actual electron densities of plasmas generated in the CTF. The fact that this technique is 
known to be the most accurate method of non-interfering analysis of plasmas 
(Huddlestone and Leonard, 1965), however, suggests that the results produced by this 
program should provide a reasonable estimate of the actual values. 

The results of the three data sets agree with what is to be expected. Variations in 
pressure on the plasma jet produce corresponding variations in particle densities and 
temperatures; this 1s particularly evident when comparing CTF146 to the other two data 
sets. Additionally, the uniform radial electron density distribution agrees with Jahn’s 
(1968) statement that an MPD cathode jet has a generally invariant radial electron 
distribution. 

A. FUTURE TOPICS 

In addition to half width measurement, actual matching of observed line profiles 
to theoretical Stark profiles is a proven method of particle density determination with an 
accuracy much greater than half width matching alone (Vidal, Cooper and Smith, 1972). 
A possible way of doing this would be by taking the Fourier transform of the line profiles 
and de-convolving them with the known Doppler profile for the plasma. The de- 
convolved profile could then be matched with known line profiles to achieve an accurate 
measurement of electron densities. This method requires accurate knowledge of the 
electron temperature, as well as background intensity levels determined in a manner other 
than presented in this thesis. 

This thesis was limited to three data sets only, which were taken from inside the 
cathode jet of an MPD test facility. Additional regions of exploration within the plasma 
jets generated in an MPD thruster could be near the surface of the anode, within the anode 
jet, in the regions between the two jets, and in the sheath region around the cathode. The 


latter instance would require a translating instrument with far finer resolution than the 


oo 


micrometer used in this experiment, being that the sheath is on the order of several 
nanometers thick. 

Additionally, a means of direct measurement of electron densities and 
temperatures in the vacuum chamber should also be explored as a way of providing 


independent verification for the analysis presented here. 
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APPENDIX A. CODE 


Following is the Interactive Data Language (IDL) code developed to determine the 
electron densities and temperatures for the data sets taken. It consists of thirteen separate 
IDL procedures, and is run by making a call in IDL to a single procedure called 
LINBROAD, called as such: 

LINBROAD, ’<filename>’, dx,fwhm,dr,Ir,dw,lw,Te,N 


Table A.1 defines the parameters as listed in the above statement: 
















dr Ratio of Doppler half width to full 
half width 
Ratio of Lorentzian half width to full 
half width 


rs 


Table A.1. Definition of Input/Output Parameters for IDL Code. 


The names of the parameters are entirely up to the user, within the restrictions of filename 


size and restricted keywords for the IDL language. 


yf 


; LINBROAD 
; Master procedure for analysis of broadened 486.133 nm hydrogen line 


; Input Variables: 


ctfno Name of output file for data storage 


;Output Variables: 


: dx Position of scan (mm) 

‘ fw full half-width of line 

dr Ratio of Doppler half width to observed half width (full) 

: lr Ratio of Lorentzian half width to observed half width (full) 
dw Doppler width, angstroms 


: lw Lorentzian (Stark) width, angstroms 
: Te Electron temperature, K 
: N Electron density, cm4-3 


pro linbroad,ctfno,dx,fw,dr,lr,dw,lw,Te,N 


fitfile,Ir,dr,fw,dw,lw,dx ;Call procedure FITFILE 
NeTe,fw,lw,dw,Te,N ;Call procedure NETE 
Storit,ctfno,dx,fw,dr,Ir,dw,lw, Te,N ;Call procedure STORIT 
retum ;Retumn 

end ;End procedure 

@fitfile ;Compile FITFILE 
@NeTe ;Compile NETE 

@storit ;Compile STORIT 
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‘ FITFILE 


; Procedure sequentially reads in ‘scan1' through ‘scan31', preps and centers data, and 
; determines the Doppler and Stark widths of each observed line profile. 


;Output Variables 


: lr Ratio of Lorentzian half width to observed (full) half width 
dr Ratio of Doppler half width to observed (full) half width 
fw Full half width of observed profile, angstroms 

dw __— Doppler width of observed profile, angstroms 

: lw Lorentzian (Stark) width of observed profile, angstroms 

: dx Position of scan (mm) 


pro fitfile,Ir,dr,fw,dw,lw,dx 


junk = check_math(0,0) 


fw = dblarr(3 1) ;Define output arrays fw, dw, lw, Ir, dr 
dw =fw 

lw = fw 

lr=lw 

r= ir 

dx = (findgen(30) - 15)*0.083*2.54 ;Generate dx 


for m = 0,30 do begin 


case 1 of ‘Determine scan to read 


(mlt15): filename = strcompress(‘scan'+string(3 1-m),/remove_all) 
else : filename = strcompress(‘scan'+string(m-14),/remove_all) 


endcase 
getfile,filename,wl,int,par, fit ;Call procedure GETFILE 
delta = sqrt(par(3)) ‘Initial Doppler width estimate (half width of 
;Lorentzian fit performed in GETFILE 
lint = int-par(0) ;Adjust scan and fit for offset determined in 
GETFILE 


lfit = lfit-par(O) 


ao 


a= fltarr(4) ;Define imnitial-guess array for Voigt 


parameters 

a(O) = lint(1500) Peak intensity 

atl =292 Damping parameter 

a(2) = -delta Doppler width 

a(3) = -par(0) Offset parameter 

b=a tore initial guesses 

w = fltarr(n_elements(wl))+1 ;Define weight vector (no weighting) 

chsq = 999 sset Chi-Squared to 999 

chsqp = 1001 ;Define previous Chi-Squared value at 1001 
vfit = dblarr(3000) ;Define array to store fit of Voigt profile 
while chsq It chsqp do begin ;Execute loop until Chi-Square minimized 

chsqp = chsq ;Update chisqp 

pre = a(2) ;Store current Doppler width 

vpre = vfit ‘Store current fit 


vfit = curvefit(wl,int,w,a,function_name = 'voigt_a',/noderivative) ;Fit Voigt profile 


vint = int + a(3) ;Adjust intensity by offset from fit 
vfit = vfit + a(3) ;Adjust fit 
dop = a(2) ;store Doppler width 
chsq = total((vint-vfit)*2/vfit) ;Compute Chi-Squared of fit 
a= Dp ;Reset a with initial values, using Doppler 
a(2) = dop ;width from last fit 
endwhile ‘-Re-execute fit 
vfit = vpre ;Best fit 
halfwidth,wl,vfit,del,ind ;Call procedure HALFWIDTH 
getpar,pre,del,damp,l,d ;Call procedure GETPAR 
Bern) al ;Lorentz width fraction 
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dr(m) =d ;Doppler width fraction 


prowidth,wl,lint,fwhm ;Call procedure PROWIDTH 
dw(m) = d*fwhm ;Doppler width 

lw(m) = 1*fwhm ;Lorentzian (Stark) width 
fw(m) = fwhm sFull half width 


print,{dw(m),lw(m),fw(m)] 

endfor 

return ‘Return 

end sEnd 

@ getfile ;Compile called procedures 
@ getdata 

@voigt_a 


@halfwidth 
@prowidth 
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; GETFILE 


- Reads appropriate scan, and fits an initial Lorentzian profile to determine line center and 
; provide initial estimate parameters for fitting a Voigt profile. 


‘Input Variables 

: filename Name of scan to be read 
;Output Variables 

lam Wavelength vector, angstroms 
: int Intensity vector, volts 

par Lorentzian fit parameters 

fit Lorentzian fitted profile 


pro getfile,filename,lam,int,par,fit 


wl 1=dblarr(5000) ;Define arrays for raw data 

intl=wl1 

wl0=wl1(0) ;Set input variables for reading data 

intO=int1(Q) 

wl=dblarr(4600) ;wl = wavelength 

val=wl ;vai = voltage 

st=" ;Blank string 

openr, |, filename ;Open input file 

for j=0,5 do readf,1,st kip first six lines 

for j=0,4999 do begin ;Read in 5000 point wavelength and voltage 
vectors 

readf, 1 ,wl0,intO 

wl1Q)=wlO 

int1()=intO 

endfor 

close, 1 ;Close input file 

wl=wl11(200:4799) ;Delete first and last 200 points 


val=int 1(200:4799) 


profcor,wl,val,lam1 ,int1,coef,fit1 ;Call procedure PROFCOR 
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low = where(lam1 It coef(2)-3.0,locount) set cutoff at +/- 3 angstroms from line center 


high = where(lam1 gt coef(2)+3.0,hicount) 
locut = locount 
hicut = n_elements(lam1 )-hicount 


lam = extrac(lam1,locut,hicut-locut) ;Extract working vectors from withing cutoffs 
int = extrac(int1,locut,hicut-locut) 
fit = extrac(fit1 ,locut,hicut-locut) 


par = coef ‘Parameters from Lorentzian fit 

dl = abs(lam(0:1499) - coef(2)) Split profile into left and right sides 
dr = lam(1499:2999) - coef(2) 

X = sort(dl) 

lint = int(x(*)) ;’Flip” left hand array 

rint = int(1499:2999) 

hint = (lint+rint)/2 ;Average left and right hand arrays 
int(0:1499) = hint(x(*)) ;’Unflip” left hand array 
int(1500:2999) = hint 

return ‘Return 

end ‘End 

@profcor ;Compile called procedures 
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; PROFCOR 


; Fit initial Lorentzian profile to data 


Input Variables 
lam_O Wavelength 
; int_O Intensity 


;Output Variables 


lam Centered wavelength 
: int Returned intensity 
coeff Fit coefficients 

yfit Fit 


pro profcor,lam_0,int_0,lam,int,coeff, yfit 


ain = [-1.0,7.0,4861.33,0.2] ;Initial parameter guess 
lorfit,lam_O,int_O,ain,yfit,a ;Call procedure LORFIT 
lam = lam_O - a(2)+4861.33 ;Set center wavelength to 4861.33 angstroms 
lorfit,lam,int_0O,ain,yfit,a ;Fit around 4861.33 angstroms 
int = int_O 

coeff =a 

return ‘Return 

end ‘End 

@lorfit ;Compile called procedures 


; LORFIT 


- Fits user-defined Lorentzian function 


‘Input Variables 

: X X vector 

y Y vector 

ain Initial parameter guess 
;Output Variables 

yfit Fitted Lorentzian 

: aout Returned parameter vector 
; sigma_a Standard deviation vector 


pro lorfit,x,y,ain,yfit,aout,sigma_a 


w = fitarr(n_elements(x))+1 ;Define weight vector (no weighting) 
aout = ain 
yfit = curvefit(x,y,w,aout,sigma_a,function_name = ‘lorentz’) ‘Perform fit 
return ‘Return 

end “End 

@lorentz ;Compile fitted function 
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; LORENTZ 


‘Defines Lorentzian function 


Input Variables 
: X X vector 
; a Parameters (constants) 


;Output Variables 


if Returned function 
pder Array of partial derivatives wrt parameters 


° 
9 


pro lorentz,x,a,f,pder 
f = a(O) + a(1) / (x - a(2))42 + a(3)) ;Define Lorentzian function 
if n_params() ge 4 then begin ;Define partial derivatives 
pder = fltarr(n_elements(x),4) 
pder(*,0) = 1.0 
pder(*,1) = 1.0/((x-a(2))42+a(3)) 
pder(*,2) = 2*a(1)*(x-a(2))/((x-a(2))42+a(3))42 
pder(*,3) = -a(1)/((x-a(2))*2+a(3))42 
endif 


end “End 
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; VOIGT_A 


;Define Voigt function for curve fitting based on the IDL “Voigt” function 


;Input Variables 

: lam Wavelength 

: a Fit parameters 

;Output Variables 

: f Returned function 

pder Partial derivative vector (not used) 


pro voigt_a,lam,a,f,pder 


u = (lam - 4861.33)/a(2) ;Dimensionless wavelength vector 
f = a(0)* voigt(a(1),u) - a(3) ; Voigt function 

return >Return 

end “End 
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sHALFWIDTH 


; Determines full half width of Voigt profile 


‘Input Variables 
‘ lam Wavelength 
V Voigt function 


;Output Variables 


hw Half width of profile, in angstroms 
a Index of half maximum point 


pro halfwidth,lam,v,hw,a 


val = abs((max(v)/2-v)/(max(v)/2)) ;Determine half max value 
a = where(val eq min(val)) ‘Locate half max 
hw = abs(lam(1499) - lam(a)) sHalf width 

return -Return 

end ;End 
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;GETPAR 


;Determines damping ratio, Lorentzian and Doppler fractional half widths given doppler 
width 
cand half width of Voigt profile. 


Input Variables 

; dw Doppler width 

: hw Half width 

Output Variables 

‘ d Damping parameter 

: Ir Lorentzian fractional half width 

: dr Doppler fractional half width 

pro getpar,dw,hw,d,Ir,dr 

damp = fltarr(29) ;Define input arrays 
Ico = damp 

dco = damp 

damp0 = damp(0) ;Define input variables 


lcoO = Ico(0) 
dcoO = dco(0) 


openr, 1,'voigtdat’ ;Open input file 


for ) = 0,28 do begin 


readf, 1,damp0,lco0,dco0 ;Read Voigt parameters 
damp(j) = damp0 
Ico(j) = IcoO 
dco(j) = dcoO 
endfor 
close, 1 ;Close input file 
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coef = poly_fit(dco,lco,2) ;Binomial fit of Doppler and Lorentz ratios 


dr = abs(dw/hw) ;Doppler ratio of voigt profile 
Ir = coef(0)+coef(1)*dr+coef(2)*dr“%2 ;Lorentz Ratio 

d = Ir/dr ;Damping parameter 

return ‘Return 

end ‘End 
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;PROWIDTH 


; Determines full half width of data. Profile width computed by averaging data values over 


a 20 


; point spread, centered about every data point except the first and last ten. Values are then 
;compared to the average of the 20 values centered around the line center. 


Input Variables 

] Wavelength 
int Intensity 
;Output Variables 

: fw Full half width 


pro prowidth,I,int,fw 
left = 1(0:1499) 

right = 1(1500:2999) 
li = int(0:1499) 

ri = int(1500:2999) 


half = total(int(1489:1509))/42 


lcomp = dblarr(1480) 


rcomp = Icomp 


Ival = lcomp 
lind = lcomp 
rval = lcomp 
rind= lcomp 
is taken (right) 


for a=0,1479 do begin 


split profile into left and right 


;Average 20 points around line center 
sand find half maximum 


;Comparison arrays - percentage 
difference 

sbetween averaged values and half 
max 


;Average values of every 20 points in data 
(left) 
Index of point around which average is taken 
(left) 


;Average values of every 20 points in data (right) 
;sIndex of point around which average 


;Average 20 points about every point and 
record 
;Index of each point averaged 


71 


b =a+1l0 

Ilval(a) = total(li(b-10:b+10))/21 
Icomp(a) = abs((Ival(a) - half)/half) 
lind(a) = b 

rval(a) = total(ri(b-10:b+10))/21 
rcomp(a) = abs((rval(a) - half)/half) 
rind(a) = b 

endfor 

lx = where(Icomp eq min(Icomp)) 
lhw = abs(left(lind(ix)) - 4861.33) 
rx = where(rcomp eq min(rcomp)) 
rhw = abs(right(rind(rx)) - 4861.33) 
fw = rhw+lhw 


return 


end 


eZ 


;Compute index closest to half max 
(left) 
;Compute half width (left) 


;Compute index closest to half max 
(right) 

;Compute half width (right) 

‘Full half width 


‘Return 


‘End 


-NETE 


;Computes electron density and temperature 


;Input variables 


: fw Full half width 

: lw Lorentzian (Stark) width (full) 

: dw Doppler width (full) 

;Output Variables 

: Te Temperature (K) 

: N Electron density (cm‘-3) 

pro NeTe,fw,lw,dw,Te,N 
halfest,fhw,logNe,T,spln ;Call procedure HALFEST 
a = n_elements(fw) set counting variable 
In = dblarr(a) slog(N) vector 
IN = In -N vector 
d = dw/2 ;Half Doppler width 
mi 1-6/e-27 ;sHydrogen mass 
k = 1.38e-23 ‘Boltzmann constant 
c= 3e8 speed of light 
lam0 = 4861.33 ;Center wavelength 


Te = m/(2*alog(2)*k)*(-c*d/lam0)*2 
logTe = alog!0(Te) 
hwest = fltarr(a,4) 


for } = 0,a-1 do begin 


for i = 0,3 do begin 


;Compute electron temperature 
slog(Te) 


;2D array of half width versus log(N) 
sto be interpolated from fhw 


;spline interpolation to determine 
log(N) vs 
slog(dlw) for T 


hwest(j,1) = spl_interp(alog10(T),fhw(*,1),spiIn(*,1), logTe()) 


endfor 


IN(j) = interpol(logNe,alog 10(hwest(j,*)),aloglO0(lwG))) | ;Linear interpolation of log(N) 
vs 
slog(lw) to determine log(N) 


endfor 
N = 104IN ;Electron density 
return ;Return 
end ;End 
@halfest ;Compile 
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SHALFEST 


;Read in log(N) vs alpha table and generate cubic spline for interpolation of electron density. 


;Output Variables 


hw Full stark widths (2D array) 

: logN Log (base 10) of electron density vector 
i Electron temperature vector 

: Spln Spline coefficients (2D array) 

pro halfest,hw,logN,T,spln 


openr, | ,'alpha’ 
logN = fltarr(4) 
T = logN 
hw = fitarr(4,4) 
spln = hw 
b = 2.0/3.0 
for 1 = 0,3 do begin 
for } = 0,3 do begin 
readf,1,TO,N,a0 
TQ) = TO*1000 
if TO eq 5 then logNGg) = N 
hw(i,j) = a0*2.5e-9*(104N)4b 
endfor ;j 
endfor 51 


for i= 0,3 do begin 


spln(*,i) = spl_init(alog10(T),hw(*,1)) 


;Open input file 


;Define vectors and arrays 


;Electron density exponent 


‘Fill T vector and half width matrix 


;Generate splines 


endfor 

Close, 1! ;Close input file 
return ‘Return 

end ‘End 


‘STORIT 


;Writes results of analysis to file 


;sInput Variables 


filename 
; dx 

fw 

dr 

Ir 

: dw 

lw 

: Te 

N 


File name to write data to 
Scan position, mm 

Full half width 

Doppler fractional width 
Lorentzian fractional width 
Doppler width 

Lorentzian (Stark) width 
Electron temperature, K 
Electron density, cm‘-3 


pro storit,filename,dx,fw,dr,lr,dw,lw,Te,N 


openw, | filename 


for j = 0,30 do begin 


;Open write file 


printf, 1,dx(Q),fwG),drG),Jrg),dwG),lwG),TeG).NG) ;Write data to file 


endfor 
close, | 
return 


end 


‘Close write file 
>-Return 


‘End 
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APPENDIX B. TABULATED RESULTS 


On the following pages are, in tabulated format, some of the results from CTF’s 
144, 145, and 146. Included are radial position, full half width, Doppler and Stark half 


widths, electron temperature, and both the log and decimal values for the electron density. 


af 








| __—_s|talf Widths, Angstroms | TT 
R, cm Full Doppler Stark T,K logNejcm*’-3 

0.89 
0.894 
0.914 
0.898 0.7869 
0.906 
0.902 0.7808 
0.898 0.7828 
0.902] 0.3043} 0.7891 
0.906 0.7808 
0.914 
0.84328 0.898 0.7807 
0.89 0.7763 
0.85 0.7494 
0.858 
enenG 0.85 
0.21082 0.842] 0.2848] (0.7361 
0.85 0.7365 
0.846 
0.842 
0.842 
1.26492 0.846] _0.2886| _0.7378 
0.85 0.7379 
0.85 
1.89738 0.846] 0.2967 __—(0.7318 
0.846 0.7405 6744 
0.83 0.7306 
0.846 0.7455 
0.834| 0.2799] 0.7307 
2.95148 0.822) 0.2777 
3.1623 0.83] 0.2761 











Table B.1. CTF 144 Results 
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| __|Half Widths, Angstroms | ae le 
R, cm Full oF Stark T,K logNe|cm*-3 
0.7027 
-2.95148 0.2659] 0.6916 14.33 
-2.74066 0.6908 
-2,52984 0.7047 
-2.31902 0.2666] (0.7178 
| -2.1082] 0.81] 0.2874] ~—_—0.6982 
-1,89738 0.7303 14.364 
-1.68656 0.7272 14.36 
1.47574 0.7297 
-1.26492 0.2899| 0.7369 
0.7259 6818 
0.826 0.7306 5942 
0.7235 14.3565 
0.7149 14.3489 
-0.21082 0.2797| (0.7174 14.3516 
ft in| 0.2742] 0.7079 
0.7147 
0.7284 
0.63246 0.7256 
0.84328 0.83 0.7265 14.3596 
1.0541 0.7136 
1.26492 0.7329 
1.47574 0.2817| 0.7339 
1.68656 0.7242 
1.89738 0.83 0.7306 
2.1082 0.7461 
2.31902 0.7195 
2.52984 0.7233 14.357 
2.74066 0.7208 
2.95148 0.7074 
0.285 0.718 6752) 14.3514] 2.2459496 








Table B.2. CTF 145 Results. 


qe 











|__|Half Widths, Angstroms ee 
R, cm Full Dp Stark T,K logNeicm’-3 

-2.74066 3,7966451 
2.52984 14393 3.8256039 
-2.31902 15145 3.9472988 
-2.1082 0.4023 1.022 13457 
| -1.89738| 1.206] 0.4188] 1.046] 14577 38141696 
-1.68656 12995 
-1.47574 1.042 13762 3.7931498 
-1,26492 14023 3.7290654 
0.4137 14224 3.8832916 


io) 


g 1382 
0.63246 1.037 1373 
1.214 1501 


14.5894}; 3.885080 
5 


14.5765| 3.771377 


14.5838] 3.835305 


WwW I~ | oO 
CO TN 











Nn 


2.7406 


15578| 14.5716] _3.7290654 
a a ee 1.039 13174] 14.5778] 3.7826835 
0.21082{ 1.21 0.4109] ~—s1.057| 14033] ~—:14.5879] 3.8716849 
| 9.42164 1.198] 0.4035) 1.048 13537| 14.5833] 3.8308928 
14199] ___14.584]_3.8370725 
0.84328 13788] 14.5788] 3.7914034 
1.0541 14199] 14.5676] 3.6948771 
1.26492] 1.182] 0.4025 1.031] ~——s-13462| ~——«:14.573] 3.7411059 
1.47574 1.174] «0.3941 3.7290654 
1.68656] —«:1.186{ 0.4029 ~—=s1.036] ~—s:13491| 14.5755] 3.7627035 
1.89738 1.19] 0.416 ——s1.03) ~=——:14379| 14.5718] 3.7307831 
2.1082 1.194 0.416) ~—s1.035| 14373] 14.5746] 3.754914 
2.31902] —-:1.198} 0.4062] ~—s«1.046) ~—s-13713] 14.5821} 3.8203223 
2.52984 1.206 0.4141) ~—si.05| «14285 3.8335 
g 1.045 


0.4083 


Table B.3. CTF 146 Results 


14.5811} 3.811535 
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APPENDIX C. DATA TABLES 


Two tables are presented here which were used in the analysis program listed in 
Appendix A. Table C.1, from Huddlestone and Leonard (1965) is a tabulation of Voigt 
profile characteristics, where A, is the Doppler half width, A, is the Lorentzian half width, 
and A,, is the observed half width of the Voigt profile. The first column, A,/Ap, 
represents the damping parameter of the profile. The second and third columns display, 
respectively, the ratios of Lorentzian and Doppler half widths to the profile width. Using 
columns two and three, a second order polynomial fit of A, vs Ap was generated, from 
which a value of A, could be determined, given Ap and i,,. 

Table C.2 from Griem (1974) lists computed «,, relative to electron densities and 
temperatures. These data are interpolated using cubic splines in order to determine the 
half widths corresponding to each electron density for the value of 7, computed from the 
Doppler half width of the profile. Once the half widths are determined, a linear 
interpolation (on a log-log scale) of NV, vs half width is made for the Lorentzian half width 


of the profile, and an electron density is determined. 
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Table C.1. Half widths of Voigt profiles (From Huddlestone and Leonard, 1965) 
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O1 


| tog(Ne) 14415 
Electron Density, cm“-3 |Full half width, alpha (unitless 
5000 7.62e-0 7.87e-02] 8.08e-02| 7.65e-02 


NO 


7.90e-02| 8.14e-02 
sid Full half width, angstroms | 
| 5000 __—_—=*4.0,4104198] _1.9675|9.376009 


© 


| 


4 
10000 0.4168831 2.007519.7473366145.835598 
20000 0.4184989 2.0375] 9.979416]48.582502 
40000 0.4255009 2.035] 9.99102]49.929024 


Table C.2. Half width of theoretical Stark profiles vs N,, 7,. (From Griem, 1974). 


- 
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